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NOTATION 


Strength  member  cross-sectional  area 
Constant  coefficient  in  hydrodynasiic  loading  function 
Coefficient  of  cos(n+)  terms  in  hydrodynamic  loading  function 
Cable  buoyancy  it:  fluid 

Coefficient  of  sin(n+)  terms  in  hydrodynamic  loading  function 
Cable  normal  drag  coefficient  based  on  frontal  area 
Cable  side  force  coefficient  based  on  frontal  area 
Cable  diameter  or  thickness 

Modulus  of  elasticity  at  any  point  on  a  stress-strain  curve 
Slope  of  stress-strain  curve  at  aero  stress 

Coefficient  of  stress  term  in  equation  for  modulus  of  elasticity 
Cable  hydrodynamic  force  component  per  unit  length  normal  to  a  cable 
element  in  the  plane  defined  by  the  cable  element  and  the  free 
stream  direction 

Cable  hydrodynamic  side  force  per  unit  length  at  the  orientation  of 
the  cable  that  produces  the  maximum  value 

Body  or  cable  force  components  in  the  X-,  Y- ,  and  Z-directions 
Pode  frictional  parameter 

Normal,  tangential,  and  side  force  hydrodynamic  loading  functions' 
Cable  hydrodynamic  force  component  per  unit  length  tangential  to  a 
cable  element 

Cable  hydrodynamic  force  component  per  unit  length  normal  to  the 
plane  defined  by  a  cable  element  and  the  free-stream  direction 


NOTATION  (Continued) 


Cable  dreg  per  unit  length  when  the  cable  it  normal  to  the  free- 

stream  direction 

Reynolds  number 

Stretched  cable  length 

Unstretched  cable  length 

Cable  length  printout  interval 

Cable  tension 

Free-stream  velocity 

Cable  weight  per  unit  length  in  a  fluid 
Cable  weight  per  unit  length  in  air 

Axis  of  space-fixed  orthogonal  coordinate  system  positive  in 
direction  of  tow  (negative  in  free-stream  direction) 

Axis  of  space-fixed  orthogonal  coordinate  system  positive  to  the 
right 

Axis  of  space-fixed  orthogonal  coordinate  system  positive  in  the 
direction  of  gravity 

Intermediate  coordinate  system  defined  by  a  rotation  of  angle  8 
about  the  X-axis 

Distances  along  the  X,  Y,  and  Z  directions 

Cable  kite  angle  measured  from  the  Z-axis  to  the  tangent  of  a  cable 
element  projected  onto  the  Y-Z  plane 
Cable  strain 

Fluid  kinematic  viscosity 
Flui'1  density 
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ABSTRACT 


A  generalised  computer  program  for  predicting 
the  three-dimensional ,  static  configurations  and 
tensions  of  extensible ,  flexible  cable  systems  is 
described.  The  program  is  written  in  BASIC  lan¬ 
guage  for  use  on  desktop  computers*  The  report 
includes  development  of  the  equations*  description 
of  the  program,  instruction  for  program  usage,  sam¬ 
ple  problems,  and  a  program  listing  for  a  Hewlett- 
Packard  9836  desktop  computer. 


ADMINISTRATIVE  INFORMATION 

The  work  described  in  this  report  was  performed  in  support  of  a  number  of 
projects  sponsored  by  the  Office  of  Naval  Research,  the  Naval  Air  Systems  Command, 
the  Naval  Sea  Systems  Command,  and  the  Naval  Electronics  Systems  Coraand.  This 
report  was  written  under  Office  of  Naval  Research  Uork  Request  30037  of  S  January 
1983,  David  Taylor  Naval  Ship  Research  and  Development  Center  Work  Unit  1541-202. 


INTRODUCTION 

This  report  describes  a  generalized  computer  program  for  predicting  the 
three-dimensional,  static  configurations  and  tensions  of  extensible,  flexible 
cable  systems  in  a  uniform  stream.  The  program  is  written  in  BASIC  language  and 
is  specifically  intended  to  provide  time-efficient  solutions  on  a  desktop  compu¬ 
ter.  Although  the  program  initially  was  written  for  use  with  dwwlett-Packard  (HP) 
computers,  it  can  be  adapted  with  suitable  variations  in  input/output  format  to 
other  computers  that  use  BASIC.  To  date,  the  program  has  been  used  on  an  HP-9826, 
an  HP-9835,  an  HP-9836,  and  an  HP-87.  A  program  listing  for  the  HP-9836  is 
contained  in  the  Appendix. 

The  program  is  designed  to  solve  initial  value  problems  when  the  conditions 
at  one  end  of  the  cable  are  known.  Boundary  value  problems  must  be  solved  by  an 
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iterative  procedure  to  eetiafy  prescribed  conditions  at  both  ends  of  the  cable* 
The  prograa  allows  the  cable  to  assume  any  configuration  without  regard  to  or 
prior  knowledge  of  cable  quadrant.  The  correct  configuration  is  automatically 
established  by  the  signs  of  the  initial  force  components.  Cable  elasticity  is 
entered  as  an  option.  In  some  cases,  as  for  example  with  steel  cables,  elongation 
may  be  insignificant.  If  this  is  the  case,  program  input  is  simplified  by 
choosing  the  nonelestic  option.  If  elongation  is  significant,  as  with  many  com¬ 
posite  end  synthetic  lines,  then  the  elastic  option  will  provide  greater 
accuracy.  To  use  this  option,  cable  elastic  modulus  and  cross-sectional  area  must 
be  entered  as  input  data.  The  elastic  modulus  may  be  entered  either  as  a  constant 
or  as  a  linear  function  of  cable  stress.  Synthetic  lines  frequently  exhibit  e 
variable  isodulus  proportional  to  stress. 

The  program  allows  intermediate  bodies  to  be  added  at  any  decired  point  along 
the  cable.  In  addition,  certain  properties  may  be  changed  along  the  cable  as 
desired;  for  example,  the  cable  stay  be  part  ribbon-faired  and  part  bars  or  the 
fluid  My  change  from  water  to  air.  Trail,  kite,  depth,  and  cable  angle  stop 
values  also  may  be  imposed.  To  include  an  intermediate  body  or  to  change  a  cabls 
parameter,  the  cable  length  to  that  point  or  a  stopping  parameter  it  specified. 
The  program  then  will  calculate  the  cable  configuration  to  the  specified  value  and 
request  the  option  desired.  In  this  way,  the  computation  can  be  continued  to  any 
desired  point  of  termination. 

This  report  includes  derivations  of  the  generalized  equations  which  govern 
the  static  configuration  and  tension.  Also,  the  assumed  representations  of  the 
fluid  forces  and  cable  elastic  properties  are  discussed  in  detail.  Finally;,  the 
computer  program,  including  the  method  of  solution  and  its  usage,  i9  described. 
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DERIVATION  OF  CABLE  CONFIGURATION  EQUATIONS 


The  differential  equations  describing  the  three-dimensional  static  configura¬ 
tion  and  tension  of  a  cable  in  a  uniform  stream  are  derived  from  the  equilibrium 
cc  external  forces  acting  on  an  element  of  the  cable.  A  free-body  diagram  showing 
a  segment  of  cable  of  elemental  Length  ds  acted  upon  by  hydrodynamic,  hydrostatic, 
gravitational,  and  tension  forces  is  illustrated  in  Fig.  1.  The  (X,  Y,  Z)  coor¬ 
dinate  system  shown  is  a  right-hand,  orthogonal  system  fixed  in  space  with  the 
X-axis  positive  in  the  direction  of  tow  (or  negative  in  the  free-stream  direction) 
and  the  Z-axis  positive  in  the  direction  of  gravity. 

The  equations  are  conveniently  derived  for  an  orthogonal  coordinate  system 
fixed  to  the  cable.  In  Fig.  1,  the  hydrodynamic  force  has  been  resolved  into 
components  Fds,  Gds,  and  Hds  where  F  is  the  force  component  per  unit  length  normal 
to  the  cable  in  the  plane  defined  by  the  cable  element  and  the  free-stream  direc¬ 
tion,  G  is  the  force  component  per  unit  length  tangential  to  the  cable,  and  H  is 
the  (side)  force  component  per  unit  length  normal  to  the  plane  defined  by  the 
cable  element  ;'.nd  the  free-stream  direction.  The  cable-fixed  coordinate  system 
defined  by  the  directions  of  F,  G,  and  H  may  be  obtained  by  first  rotating  the 
spatial  system  by  an  angle  8  about  the  X-axis  and  then  rotating  the  resulting 
intermediate  (Xp  Yp  Z^)  coordinate  system  by  an  angle  4>  about  the  Ypaxis.  In 
towed  systems  nomenclature,  the  angles  8  and  $  generally  are  referred  to  as  the 
kite  angle  and  the  cable  angle,  respectively.  Since  the  orientation  of  the  cable 
changes  in  space,  the  angles  8  and  $  are  functions  of  cable  length  s. 

For  equilibrium,  the  sums  of  the  force  components  in  each  of  the  directions 
F,  G,  and  H  must  all  be  zero.  A^  shown  in  Fig.  1,  the  sum  of  the  force  components 
in  the  direction  of  F  is 
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(1) 


Fds  +  (Wft  -  B)  cosfl  cos*  ds 
-(T  +  ijdT)  sinC*  +  Ijd*)  costedS)  cos* 

+(T  ♦  IjdT)  cos(*  +  ljd*)  sin* 

+(T  -  *5dT)  sin(*  -  %d*)  cos  (tedB)  cos* 

-(T  -  4dT)  cos(*  -  ^d4> )  sin*  =  0 

where 

F  is  the  hydrodynamic  force  component  per  unit  length  normal  to  a  cable 

element  in  the  plane  defined  by  the  cable  element  and  the  free-stream 
direction, 

s  is  the  stretched  cabie  length* 

W  is  the  cable  weight  per  unit  length  in  air, 

B  is  the  cable  buoyancy  per  unit  length, 

T  is  the  cable  tension, 

B  is  the  kite  angle  measured  from  the  Z-axis  to  the  tangent  of  the  cable- 
element  projection  onto  the  Y-Z  plane,  and 
*  is  the  cable  angle  measured  from  the  free-stream  direction  (X-axis)  to 


the  tangent  of  the  cable  element. 

From  trigonometry, 

sin(*  +  4d<J>)  =  sin*  costed*)  +  cos*  sinted*)  (2a) 

cos(*  +  4d*)  =  cos*  costed*)  -  sin*  sinted*)  (2b) 

sin(*  -  *sd*)  =  sin*  costed*)  -  cos*  sinted*)  (2c)  j 

I 

cos(*  -  ^sd*)  =  cos*  costed*)  +  sin*  sinted*).  (2d) 

Also,  s'nce  d*  and  dB  approach  zero, 

costed*)  =  1,  (3a) 

sinted*)  =  ^d*,  (3b) 
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gives 


coa(^dS)  *  1,  and  (3c) 
cos(-%d8)  *  1.  (3d) 
Substituting  Eqs.  2  and  3  into  Eq.  1  and  then  expanding  and  combining  terms 


Fds  +  (Wa  -  B)  cosB  cos*ds  ~  Td*  *  G 
or 

-T  4^  ♦  F  ♦  (W  -  B)  cosS  cos*  *  0. 
ds  a 

The  sum  of  the  force  components  in  the  direction  of  G  is 

Gds  -  (Wa  -  B)  cos0  sin*ds 

+(T  ♦  IjdT)  sin(*  +  ^d*)  cos^dB)  sin* 

♦  (T  +  *jdT)  cos(*  +  ^d*)  cos*  (5) 

-(T  -  JjdT)  sin(*  -  *id*)  cos  (-*id$)  sin* 

-(T  -  ^dT)  cos(*  -  *Sd*)  cos*  *  0 

where  G  is  the  hydrodynamic  force  component  per  unit  length  tangential  to  the 
cable  element. 

Substituting  Eqs.  2  and  3  into  Eq.  5  and  then  expanding  and  combining  terms 

gives 

Gds  -  (W„  -  B)  cosB  sin*ds  +  dT  *  0 
a 

or 

4~  +  G  -  (W  -  B)  cosB  sin*  =  0. 

ds  a 

(6) 

The  sum  of  the  force  components  in  the  direction  of  H  is 
Hds  -  (Wg  -  B)  sinB  ds 

-(T  +  JjdT)  sin(*  +  4d*)  sin^dB)  (7) 

+(T  -  ^dT)  sin(*  -  %d*)  sin(-iidB)  *  0 
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where  H  is  the  hydrodynamic  force  component  per  unit  length  normal  to  the  plane 
defined  by  the  cable  element  and  the  free-stream  direction. 


Since  dfl  approaches  zero, 

sin(*sdB)  *  *idfl  and  (8a) 

sinC-^dfl)  =  -*sd8.  (8b) 

Substituting  Eqs.  2,  3  and  8  into  Eq.  /  then  expanding  and  combining  terms 


gives 

Hds  -  (Wfl  -  B)  sinBds  -  Tsin*d8  *  *«dTdBd*cos*  =  0 
or,  since  the  third  order  term  is  negligible, 

-T  sin*  4^  +  H  -  (W  -  B)  sinfl  =  0. 
as  a 


(9) 


Equations  4,  6,  and  9  are  the  three  differential  equations  of  force 

equilibrium  in  the  cable-fixed  reference  system.  The  differential  equations  for 
cable  displacements  in  the  X,  Y,  and  Z  directions  can  be  obtained  from  Fig.  1  by 
inspection.  They  are 

dx  =  ds  cos* 

dy  =  -ds  sin<t>  sinfl 
dz  =*  -ds  sin*  cosfl 

where  x,  y,  and  z  are  distances  along  the  X,  Y,  and  Z  directions,  respectively. 

FLUID  HYDRODYNAMIC  FORCE  COMPONENTS 

The  fluid  hydrodynamic  force  components  Fds,  Gds,  and  Hds  arise,  in  the  cases 
of  F  and  G,  from  fluid  drag  forces  and,  in  the  case  of  H,  from  fluid  side  or  lift 
forces. 


or 


or 


or 


~  =  cos*, 

(10) 

4^  =  -sin* 
ds 

sinfl,  and 

(11) 

dz  .  , 

^  *  -mt 

cosfl 

(12) 
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FLUID  DRAG  FORCES 


The  fluid  dreg  force  components  per  unit  cable  length  F  and  G  lie  in  the 
plane  defined  by  a  cable  element  and  the  free-stream  flow  direction.  Generally 
these  components  are  assumed  to  be  the  product  of  a  drag  per  unit  length  that  is 
a  function  only  of  Reynolds  number  and  a  loading  function  that  is  a  function  cnly 
of  cable  angle  Under  this  convention,  the  normal  force  per  unit  length  F  and 
the  tangential  force  per  unit  length  G  are  of  the  following  form: 

F(Re,*)  =  -R(Re)  •  fnU)  (13) 

G(Re,4>)  =  -R(Re)  •  ftU)  (14) 

where 

R  is  the  cable  drag  per  unit  length  when  the  cable  is  normal  to  the 

free-stream  direction  (R  =  *jpCRV2d), 

Re  is  the  Reynolds  number  (Re  *  Vd/v), 
p  is  the  fluid  density, 

CR  is  the  cable  normal  drag  coefficient  based  on  frontal  area  expressed 
as  a  function  of  Reynolds  number, 
d  is  the  stretched  cable  diameter, 

V  is  the  free-stream  flow  velocity, 

v  is  the  fluid  kinematic  viscosity,  and 

fn,ft  are  the  normal  and  tangential  hydrodynamic  loading  functions, 
respectively. 

The  negative  (-)  signs  in  Eqs.  13  and  14  appear  since  the  cable  drag  is  in 
the  negative  (-)  X  direction  as  defined  in  Fig.  1. 


The  hydrodynamic  loading  functions  fn  t  are  commonly  expressed  in  the  form  of 
a  trigonometric  series*’*  a.*  follows: 

00 

f  ($)  *  A  +  A  cos(n$)  ♦  B  sin(n$)  (15) 

n,t  °  "l  n  n 

where  A0 ,  An,  and  BR  are  constant  coefficients.  Generally,  cable  loading  is 
adequately  expressed  by  inclusion  of  only  the  first  five  terms  (n  *  2)  of  Eq.  15. 

FLUID  SIDE  FORCES 

Fluid  side  forces  occur  as  a  result  of  cable  instabilities  or  lateral  shape 
asymmetries.  In  the  case  of  twisted  wire  ropes  or  double-armored  electro¬ 
mechanical  cables,  the  helical  twist  of  the  exposed  wire  strands  causes  a  lateral 
shape  asymmetry  which  can  produce  significant  side  forces. 

In  a  manner  similar  to  that  used  with  drag  forces,  the  side  force  component 
per  unit  length  H  is  expressed  conveniently  as  the  product  of  a  side  force  that  is 
a  function  only  of  Reynolds  number  and  a  loading  function  that  is  a  function  only 
of  cable  angle  $  as  follows: 

H(Re,  <p)  =  Fg(Re)  •  fg(<|>)  (16) 

where 


Fg  is  the  cable  side  force  per  unit  length  at  the  orientation  of  the  cable 
that  produces  the  maximum  value  (F  =  4pC_V*d), 

S  9 


Cg  is  the  cable  side  force  coefficient  based  on  frontal  area  expressed  as  a 
function  ot  Reynolds  number,  and 
fg  is  the  hydrodynamic  side  force  loading  function. 
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The  maximum  side  force  Fs  may  or  may  not  occur  when  the  cable  is  normal  to 
the  f ree— stream.  The  side  force  loading  function  f^  also  is  adequately  expressed 
using  the  first  five  terms  of  the  trigonometric  series  indicated  in  Eq.  IS. 


CABLE  ELONGATION 

The  cables  used  in  towing  and  mooring  applications  generally  exhibit  three 
types  of  elongation  as  follows: 

1.  Permanent  Elongation  —  This  type  is  nonrecoverable  and  is  associated 

primarily  with  fabrication  looseness  inherent  in  new  cables,  ropes,  and 
lines  before  the  wires  or  strands  are  seated.  The  amount  of  permanent 
elongation  depends  on  the  construction  type  and  the  degree  of  pre¬ 
stressing.  Permanent  elongation  becomes  fairly  constant  after  repeated 
loadings.  For  wire  rope,  this  type  of  elongation  is  typically 

between  0.25  and  0.33  percent  of  initial  length  depending  upon  maxia  m 

3 

stress . 

2.  Creep  Elongation  —  Th’’  s  type  of  elongation  is  associated  with  material 
properties  and  depends  on  both  the  tensile  stress  and  the  length  of  time 
the  stress  is  applied.  Creep  elongation  is  recoverable  over  a  period  of 
time  after  the  stress  is  reduced  or  removed.  Metallic  cables  and  certain 
forms  of  composite  lines  exhibit  minimal  creep  within  normal  working 
stresses,  whereas  many  synthetic  lines  exhibit  noticeable  creep. 

3.  Elastic  Elongation  —  This  type  of  elongation  is  immediately  recoverable 
and  depends  on  the  type  of  construction  as  well  as  the  material  proper¬ 
ties.  Elastic  elongation  is  essentially  proportional  to  stress  in  most 
metallic  cables  and  ropes  but  is  frequently  a  nonlinear  function  of 


I 

I 

I 

I 


I 

I 


I 


I 

I 

I 


stress  in  synthetic  lines  -  most  notably  nylon,  polyester,  end 
polypropylene. 

Of  the  three  types  of  elongation,  only  elastic  elongation  is  examined  in  this 
report  and  included  in  the  computer  program.  Permanent  elongation  can  be 
accounted  for  by  assuming  e  preconditioned  cable.  Creep  elongation  is  difficult 
to  treat  because  a  stress-time  history  is  required. 

Elastic  elongation  curves  for  representative  towline  types  are  shown  in 
Fig.  2.  These  curves  are  based  on  vendor  data  and  are  intended  to  illustrate 
general  relationships  rather  than  to  provide  specific  design  information.  Ideally 
the  elastic  properties  of  each  towcable  should  be  verified  individually.  From 
standard  mechanics  of  deformable  media^,  stress  is  related  to  elastic  strain  by 
the  equation 


where 


(17) 


Et  is  the  slope  of  the  stress-strain  curve  at  any  point  and  generally  is 
referred  to  as  the  modulus  of  elasticity, 
a  is  the  cable  stress  (o  =  T/Ac), 

T  is  the  cable  tension, 

Ac  is  the  cable  or  strength-member  cross-sectional  area, 

s-s 

e  is  the  cable  strain  (e  =  - ), 

s 

o 

s  is  the  stretched  cable  length,  and 
sQ  is  the  unstretched  cable  length. 

The  cable  modulus  of  elasticity  in  Eq.  17  depends  on  cable  construction  and, 
therefore,  should  not  be  confused  with  material  modulus  of  elasticity  which 
generally  is  higher.  For  the  nonlinear  stress-strain  curves  of  Fig.  2,  the 
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Modulus  of  elasticity  Efc  can  b«  approximated  by  a  linear  function  of  atress  o  as 
follows! 

Bt  »  EQ  ♦  Ej  o  (18) 

where 

EQ  is  the  slope  of  the  stress-strain  curve  at  a  stress  o  equal  to  aero,  and 
Ej  is  a  proportionality  faccor. 

Substituting  Eq.  18  into  Eq.  17  and  solving  for  strain,  the  result  is, 


de  = 


do 

(E  +  E,  o) 
0  1 


c 


dT 

(En  ♦  E.  T/A  ) 
o  I  c 


(19) 


Equation  19  can  be  integrated  between  aero  and  any  value  of  T  to  obtain  the 
strain  or  stretch  at  that  value.  The  result  is"* 

«  -  s_  (20a) 


K_ 

A 


or 


3  *  s  (1  ♦  £-) 
0  Ac 


where 


or 


K 


=  0 


(20b) 


(20c) 


Equation  20a, c  is  the  classical  linear  stress- strain  relationship.  For  the 
nonlinear  stress-strain  curves  suevn  in  Fxg.  2,  Eq.  20a, b  predicted  the  stretched 
cable  length  to  within  ±0.5  percent  at  all  stress  levels. 

In  terms  of  an  incremental  cable  segment  of  unstretched  length  dsQ,  Eq.  20a 
becomes 

ds  =  ds  (1  +  7-)  (21) 

o  A 


12 


STRAIN,  IS  SqI/So 

Fig.  2.  Elastic  elongation  of  various  towline  types. 
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A*  the  ceble  stretches  under  tension,  both  its  diameter  and  its  stretched 
weight  per  unit  *  ;ngth  decreases.  For  a  homogeneous  and  isotropic  material, 
transverse  strain  is  reLated  to  axial  strain  by  Poisson's  ratio. ^  However,  most 
cables  are  neither  homogeneous  nor  isotropic,  and  therefore  the  assumption  is  made 
here  that  the  stretched  cable  preserve*  its  original  volume.  Under  this  assump¬ 
tion,  the  stretched  diameter  d  and  stretched  weight  in  the  fluid  per  unit 
length  W  are  related  to  their  unstretched  states  dc  and  WQ,  respectively,  by  the 
following  relationships: 


and , 


W  ■  W 

o 


1  I* 

Tn7n 

c 


TTT7rr] 


(22) 


(23) 


where  K  is  defined  in  Eq.  20b, c. 

Substitution  of  Eqs.  21,  22,  and  23  into  the  configuration  equations  allows 
solutions  of  elastically  elongated  cables  to  be  obtained  in  terms  of  unstretched 
reference  parameters. 


PROGRAM  DESCRIPTION 

The  computer  program  consists  of  a  main  input  section  and  six  subroutines. 
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MAIM  INPUT  SECTION 


Th«  input  section  provides  e  format  for  entering  fluid  density,  towing  speed 
(or  flow  speed),  initial  body  force  components  in  the  three  principal 
directions  X,  Y,  end  Z,  cable  characteristics,  cable  length  printout  interval,  and 
various  spatial  conditions  at  which  the  integration  is  to  stop.  The  initial  cable 
angle  +  and  kite  angle  6  are  automatically  determined  within  the  program  by  the 
magnitudes  and  directions  of  the  body  force  components.  This  method  of  entering 
body  force  allows  the  integration  to  proceed  without  regard  to  or  prior  knowledge 
of  the  cable  angle  quadrant.  The  program  also  will  compute  cable  configurations, 
with  a  small  starting  error,  when  the  initial  input  force  is  zero.  This  feature 
is  useful  in  towed  array  problems  where  usually  no  body  is  attached  to  the 
outboard  end  of  the  array. 

SUBROUTINE  TANGLE 

Subroutine  "Tangle"  computes  initial  tension  T,  cable  angle  *,  and  kite 
angle  B  using  the  input  force  components.  Starting  values  for  these  variables  are 
required  to  initiate  the  program  solution. 

SUBROUTINE  CRITANG 

Subroutine  "Crito.ng"  computes  the  cable  angle  $  and  kite  angle  8  when  the 
starting  input  force  is  zero.  Although  this  subroutine  is  not  required  for  the 
program  to  obtain  a  solution,  it  does  improve  the  accuracy  at  the  outboard  end  of 
the  cable  and  does  reduce  computation  time  considerably  when  the  starting  force  is 
zero. 
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SUBROUTINE  INTEGRATOR 


Subroutine  "integrator"  solves  the  differential  equations  of  force  equilib¬ 
rium,  Eqs,  4,  6,  end  9,  using  a  self-starting,  Modified  Euler's  method.  A  large 
integration  step  4s  is  used  in  the  program  to  Minimise  solution  time.  Initially 
the  step  size  is  set  equal  to  the  specified  printout  intervel  Sine  or  to  a  length 
increment  of  SO  ft,  whichever  is  less.  At  each  integration  interval,  the 
calculation  is  iterated  about  the  changes  in  cable  angle  44  and  kite  angle  40 
until  succeeding  iteration  values  of  44  and  40  are  within  10.0002  rad 
(±0.01  deg).  If,  however,  the  computation  produces  angle  chenges  greater  than 
0.09  rad  (5  deg)  or  the  computation  does  not  converge  within  10  iteration  loops, 
then  the  ‘itep  size  is  reduced,  and  the  iteration  is  repeated  for  the  smaller 
value.  By  this  technique,  the  integration  step  size  is  efficiently  adjusted  to 
provide  a  good  degree  of  accuracy  and  a  short  computation  time.  Once  the  vari¬ 
ables  T,  4,  and  0  are  determined,  the  cable  displacement  variables  x,  y,  and  z 
(Eqs.  10,  11,  and  12)  are  computed  using  the  average  calculated  values  of  4  and  0 
between  successive  integration  steps.  A  flow  chart  of  the  integration  process  is 
shown  in  Fig  3. 

1 

In  general,  the  trigonometric  series  used  to  specify  che  cable  hydrodynamic 

i 

j 

loading  functions  is  composed  of  both  sin  and  cos  terms  and  provides  the  correct 

i 

loading  value  only  in  the  cable  angle  range  4  from  zero  to  90  deg.  Therefore,  if 
a  computation  falls  outside  this  range,  a  corrected  value  must  be  used.  This 
I  correction  is  made  within  the  program  by  assuming  that  the  normal  loading  function 

5  for  all  types  of  cables  and  the  side  ^.ading  function  for  faired  cables  are 

jj  symmetrical  about  4  =  +90  deg  and  complementary  about  4  »  0  deg,  that  the 

I  tangential  loading  function  for  all  types  of  cables  is  symmetrical  about  4=0  deg 

and  complementary  about  4  =  ±90  deg,  and  that  the  side  loading  function  for 
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stranded  (twisted)  cables  is  complementary  about  both  6  =  0  dog  and  6  =  ± 
90  deg.  These  relationships  are  summarized  in  Table  1. 

The  integration  is  terminated  at  the  maximum  specified  value  of  cable  length 
Smax  or  at  a  specified  9top  value  corresponding  to  trail  distance  x,  kite  dis¬ 
tance  y,  depth  z,  or  cable  angle  4>.  After  each  integration  step,  the  program  com¬ 
pares  the  computed  value  to  the  specified  stop  value.  If  a  specified  stop  value 
has  been  exceeded,  the  program  performs  a  simple  iteration  to  converge  on  the 
specified  value.  The  iteration  works  by  halving  the  cable  length  segment  which 
brackets  the  point  until  the  desired  stop  value  is  approached  to  within  a  certain 
limiting  criterion.  The  convergence  criterion  is  0.01  ft  for  distance  variables 
and  0.0005  rad  (0.03  deg)  for  the  angular  variable. 


Table  1  -  Relationships  of  the  hydrodynamic  loading  functions 
in  various  quadrants  of  cable  angle  <S. 


loading 

Function 

Cable  Angle  Range  (deg) 

0 

-180  TO  -90 

-90  TO  0 

0  TO  +90 

+  90  TO  +180 

<t>i  - 

180  minus  101 

l<t>l 

0 

180  minus  0 

Norma! 

fn  m  - 

-fn  (Oil 

-fn  MM) 

+  fn  (01) 

+»n  M»i) 

Tangential 

ft  M>)  - 

-ft  10 1) 

+ft  MM> 

+  ft  (01 ) 

-f,  (0i) 

Side  (stranded) 

f*  w»  - 

+*s  M»ll 

-fs 

+f8  M»l) 

~h  (t>l) 

Side  (faired) 

'» w  ■ 

-<s  M>1) 

-f«  (cb,) 

+  fs  (0!) 

+f8  MM» 

SUBROUTINES  PODE,  RI BbON ,  AND  OTHER  LOADING 

Subroutines  "Pode"  and  "Ribbon”  pruvide  hydrodynamic  loading  functions  for 
bare  and  ribbon-faired  cables,  respectively.  These  subroutines  are  included  for 
convenience  since  these  types  of  cables  are  used  frequently.  Subroutine  "Other 
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loading"  allows  the  user  to  specify  the  coefficients  of  the  trigonometric  series 
corresponding  to  any  desired  loading.  Pode  and  ribbon-faired  loading  are  dis¬ 
cussed  in  the  next  paragraphs. 

Pode  Normal  and  Tangential  Hydrodynamic  Loading 

Pode  hydrodynamic  loading7  is  commonly  used  to  predict  the  performance  of 
unfaired  (bare)  cable  systems.  In  terms  of  the  previously  discussed  trigonometric 
series  (refer  to  section  on  FLUID  DRAG  FORCES),  the  normal  and  tangential  loading 
functions  fR  and  ft  are,  respectively* 
fn  =  s in^<J>  =  0.5  -  0.5  cos(2*) 

ft  =  f  (25) 

where 

f  is  the  Pode  frictional  parameter. 

Ribbon  Normal  and  Tangential  Hydrodynamic  Loading 

Normal  and  tangential  hydrodynamic  loading  functions  have  been  obtained  from 
at-sea  data  for  a  0.78-in.  (19.8-mm)  diameter  double-armored  electro-mechanical 
cable.®  The  ribbon  used  was  0.015-in.  (0.381-mm)-thick,  0.78-in.  ( 19 .8-mm)-wide 
polyurethane.  Ribbon  spacing  was  0.78  in.  (19.8  mm)  between  centers.  The  ribbon 
length  (equal  to  one-half  the  loop  length)  was  six  times  the  cable  diameter  or 
4.68  in.  (237.7  mm).  The  normal  and  tangential  loading  functions  fn  and  ft 
obtained  for  this  ribbon  faired  cable  are,  respectively, 

fn  -  0.4986  -  0.2499  cos*  +  0.2527  sin*  -  0.2487  cos  (2*)  (26) 

ft  =  -0.2255  +  0.3417  cos*  +  0.2255  sin*  -  0.0811  sin  (2*).  (27) 
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(28) 


These  loading  functions  correspond  to  a  normal  drag  coefficient, 

CR  =  5. 7467  -  0.93  log1Q  Re;  50,000  <  Re  <  120.000 

where  Re  is  Reynolds  number  based  on  diameter.  Hydrodynamic  loading  may  vary  with 
different  ribbon  geometries,  and  therefore  the  above  equations  should  be  used  with 
caution. 


Bare  and  Ribbon-Faired  Side  Force  Hydrodynamic  Loading 

i 

The  side  force  on  1  x  19  and  6  x  19  twisted  wire  ropes  has  been  measured  m  a 
wind  tunnel.^’*0  The  data,  reduced  to  side  loading  function  form,  are  shown  in 
Fig.  4.  These  data  represent  average  values  over  a  Reynolds  number  range  from 
32,600  to  43,400  for  the  1  x  19  rope  and  from  16,400  to  83,700  for  the  6  x  19 
rope. 


Fig.  4.  Hydrodynamic  side  loading  function  for  two  twisted  wire  ropes. 
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The  side  force  coefficients  C  ,  averaged  over  the  same  Reynolds  number 
ranges,  are  0.135  and  0.232  for  the  1  x  19  and  the  6  x  19  ropes,  respectively. 
Average  values  are  given  here  since  the  data  were  scattered  to  the  extent  that  a 
reasonable  Reynolds  number  dependency  could  not  be  established.  A  least- 
squares  fit  to  the  side  loading  data  resulted  in  the  following  trigonometric 
series  for  the  side  force  loading  function  fg: 

fc  =  1,4792  +  0.0841  cos*  -  3.0425  sin*  -  1.5633  cos  (2*) 

+  1.5213  sin  (2*).  (29) 
Although  this  loading  function  was  determined  on  the  basis  of  only  two  cable 
constructions,  the  general  form  probably  is  consistent  with  other  stranded  cables 
as  well,  and  therefore  it  has  been  entered  into  the  program  as  a  general  loading 
function  for  both  bare  and  ribbon-faired  cables.  The  side  force  coefficient  C  , 
however,  is  strongly  influenced  by  the  number  of  wire  strands  exposed  to  the  flow, 
by  the  twist  or  helical  angle,  and  possibly  by  the  ribbon  geometry. 

STORAGE  AND  TIME  REQUIREMENTS 

The  fully  working  program  without  comment  statements  requires  approximately 
10,000  bytes  of  storage  on  an  HP-9836}  with  comment  statements,  approximately 
13,000  bytes  of  storage  are  required.  The  length  of  time  required  to  compute  a 
cable  configuration  depends  primarily  on  the  change  in  cable  angle  and  the  ten¬ 
sion.  On  an  average,  an  HP-9836  requires  approximately  5  to  10  sec.  However,  for 
cases  in  which  the  initial  input  tension  is  near  zero,  computing  time  may  increase 
to  20  sec  or  more.  Computation  time  also  may  be  greater  on  other  computers.  For 
example,  an  HP-9835  is  approximately  three  times  slower  than  an  HP-9836}  an  HP-87 
is  approximately  ten  times  slower. 


ACCURACY 


The  large  integration  step  site  and  relatively  loose  convergence  criteria 
used  in  the  program  do  not  appear  to  significantly  affect  accuracy.  Comparisons 
made  with  other  cable  programs  »  indicate  differences  of  less  then  0.05  percent 
in  the  output  variables  for  the  same  input  conditions. 


PROGRAM  USAGE 

Data  are  entered  into  the  program  by  a  series  of  command  statements  which 
appear  on  the  CRT  screen.  The  command  statements  are  listed  in  Table  2.  The  pro¬ 
gram  prints  the  data  variables  immediately  after  they  are  entered  to  maintain  a 
running  check  of  the  values. 

COMMENTS  CONCERNING  DATA  INPUT 

The  following  comments  refer  to  command  statements  listed  in  Table  2. 

1.  Command  statement  number  1  selects  the  printout  medium. 

2.  The  title  statement  number  2  is  dimensioned  to  allow  a  title  length  of  up 
to  80  characters. 

3.  Statement  numbers  5  through  9  request  the  type  of  hydrodynamic  loading 
desired.  If  "Pode"  loading  is  chosen,  the  statement  PODE  LOADING  will  be 
printed  and  statement  number  6,  requesting  the  Pode  frictional  parameter, 
will  appear  next  on  the  CRT.  After  the  frictional  parameter  is  entered, 


the  program  will 

skip 

to  command 

statement 

number 

10.  Similarly 

if 

"Ribbon" 

loading 

is 

chosen,  the 

statement 

RIBBON 

LOADING  will 

be 

printed. 

If  neither 

"Pode"  nor 

"Ribbon" 

loading 

is  entered, 

the 

coefficients  of  a  generalized  loading  function  will  be  requested  (command 
statements  7,  8,  and  9).  Five  values  for  each  loading  function  fn, 
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Table  2.  Program  command  statements 
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£t>  and  fg  are  required.  In  each  cate*  the  values  to  be  entered 
represent  the  coefficients  of  the  following  trigonometric  series: 

fn,t,s  *  Ao  +  A1  cos ♦  +  B1  9*n4  +  A2  cos(2$)  ♦  Bj  sin(2<fr). 

After  the  coefficients  of  all  three  loading  functions  have  been  entered, 
the  functions  with  their  assigned  coefficient  values  will  be  printed  in 
the  output. 

4.  Command  statement  number  14  provides  the  option  of  selecting  either  a 
nonelastic  cable  or  an  elastic  cable.  If  the  nonelastic  option  is 
chosen,  the  program  skips  the  command  statements  (numbers  15  and  16) 
whii.i  relate  to  the  elastic  properties  of  the  cable.  If  the  elastic 
option  is  chosen,  the  ensuing  command  statements  (numbers  15  and  16) 
require  input  of  cable  (not  material)  modulus  of  elasticity  and  strength- 
member  cr  is-sectional  area.  Two  values  representing  the  coefficients  of 
the  equation 

«20  +  E1  (T/Ac), 

where  \c  is  cable  stress,  must  be  input  for  the  modulus  of 
elasticity.  For  a  linear  stress-strain  curve,  the  first  coefficient  E 

o 

is  the  slo,  j  of  the  curve  and  the  second  coefficient  is  set  equal  to 
zero.  Fo-  a  nonlinear  stress-strain  curve,  the  coefficient  EQ  represents 
the  slope  of  the  stress-strain  curve  at  zero  stress  and  the  coefficient 
Ej  represents  a  factor  proportional  to  stress.  The  sign  of  E^  must  be 
positive.  This  method  of  representing  modulus  of  elasticity  is  useful 
for  synthetic  lines  which  frequently  exhibit  nonlinear  stress-strain 
relationships.  Modulus  of  elasticity  is  discussed  further  in  the  section 
entitled  CABLE  ELONGATION. 
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Body  force  it  entered  into  the  program  (command  statement  number  17)  in 
the  form  of  components  in  the  three  principal  directions  X,  Y,  and  Z 
where  Fx  is  positive  in  the  direction  of  tow  (negative  in  the  free-stream 
direction),  F  is  positive  to  starboard,  and  F^  is  positive  in  the 
direction  of  gravity.  Entering  body  forces  in  this  manner  avoids 
confusion  with  respect  to  cable  angle  values  and  quadrants.  Body  drag  Fx 
should  be  entered  as  a  negative  (-)  number;  a  positive  (+)  value 
indicates  thrust.  The  initial  input  force  can  be  either  at  the  fixed  end 
of  the  cable  (which  implies  a  positive  value  for  Fx)  or  at  the  towed  body 
or  free  end  of  the  cable.  Care  should  be  exercised,  however,  in  starting 
a  computation  at  a  fixed  end;  in  this  case,  if  the  starting  cable  angle 
is  near  the  critical  angle,  the  resulting  cable  shape  is  exceptionally 
sensitive  to  the  force  component  values. 

If  external  body  force  components  are  input  into  the  program  as  the  ini¬ 
tial  condition,  specified  cable  length  Smax  (command  statement  number  18) 
should  always  be  positive  (+).  If  a  negative  cable  length  is  specified, 
the  initial  input  forces  represent  internal  cable  forces  rather  than 
external  forces.  The  cable  length  printout  interval  Sine  (command 
statement  number  19)  can  be  specified  as  either  positive  or  negative;  the 
proper  sign  is  selected  by  the  program. 

Four  stopping  conditions  other  chan  Smax  are  included  in  the  program 
(command  statement  numbers  20  through  23).  The  conditions  Xotop,  Ystop, 
and  Zstop  will  terminate  the  calculation  whenever  trail  (X-direction) , 
kite  (Y-direction) ,  or  depth  (Z-direction)  reach  the  specified  value  as 
measured  from  the  starting  point  at  the  initial  body.  The  condition 
Phistop  will  terminate  the  computation  when  the  specified  cable  angle  $ 


value  is  reached.  If  no  stopping  condition  it  specified,  the 
calculations  will  continue  until  Smax  is  reached* 

8.  Once  the  computation  is  completed  to  the  specified  stopping  point,  com¬ 

mand  statement  number  2A  will  appear  on  the  CRT  screen.  If  the  option 
ADD  BODY  is  chosen,  command  statement  numbers  17  through  23  will  reappear 
on  the  CRT  in  sequence.  If  the  option  ADD  CABLE  is  selected,  command 
statement  numbers  4  through  16  and  18  through  23  will  appear  on  the 

CRT.  The  ADD  CABLE  option  allows  the  cable  properties  as  well  as  the 
fluid  density  (as,  for  example,  from  water-to-air)  to  change.  If  the 

option  NEW  CASE  is  chosen,  the  program  will  print  the  final  cable  force 
components  and  return  to  command  statement  number  1  to  begin  a  new  cable 
configuration. 

9.  When  an  HP-9836,  an  HP-9835,  or  an  HP-9826  is  used,  most  of  the  input 

values  will  be  retained  unless  they  are  specifically  changed  or  the  com¬ 
puter  program  is  stopped.  The  exceptions  are  the  stopping  parameters 

(command  statement  numbers  20  through  23),  These  values  must  be  re¬ 

entered  each  time  the  command  statement  appears  on  the  CRT.  With  an 
HP-87,  however,  all  values  must  be  re-entered  each  time  a  command  state¬ 
ment  appears  on  the  CRT, 

SAMPLE  PROBLEMS 

The  following  three  sample  problems  illustrate  program  usage: 

Problem  1  -  Submarine  Towing  a  Submerged  Buoy 

Determine  the  length  of  double-armored  steel  cable  necessary  to  tow  a  subsur¬ 
face  buoy  300  ft  above  the  deck  of  a  submarine  at  10  kn.  Ribbon  fairing  is 


attached  to  the  cable  for  the  first  300  ft  of  length  starting  at  the  buoy  to 
reduce  cable  strumming.  Printout  of  cable  variables  is  desired  at  50-ft  intervals 
along  the  cable  length.  The  cable  will  be  assumed  to  be  nonelastic.  The  charac¬ 
teristics  of  the  fluid,  cable,  and  buoy  are  listed  in  Table  3;  program  output  is 
listed  in  Table  4, 

Problem  2  -  Ship  Towing  a  Paravane 

Determine  the  length  of  faired  towline  necessary  to  tow  a  paravane  500  ft  to 
the  side  of  a  ship  at  25  kn.  The  towline  strength  member  is  parallel-strand 
fiberglass  in  a  epoxy  matrix.  The  ship  towpoint  is  15  ft  above  the  water  sur¬ 
face.  Printout  of  cable  variable:;  is  desired  at  50-ft  intervals  along  the  cable 
length.  The  characteristics  of  the  fluid,  towline,  and  paravane  are  listed  in 
Table  5;  final  program  output  is  listed  in  Table  6.  Note  that,  since  the  towline 
passes  through  the  air-water  interface  but  the  total  tide  distance  to  the  ship  is 
specified,  the  problem  is  boundary  valued  and  must  be  iterated  to  find  the  side 
distance  at  which  the  towline  pierces  the  water  surface. 

Problem  3  -  Helicopter  Towing  Two  Depressors 

A  helicopter  is  required  to  tow  two  depressors  spaced  50-ft  apart  along  a 
braided  nylon  towline  at  a  speed  of  20  kn.  Determine  the  length  of  towline 
necessary  to  place  the  lower  depressor  at  a  depth  of  75  ft  if  the  helicopter  is 
flying  at  an  altitude  of  100  ft.  The  towline  is  unf aired.  Both  depressors 
produce  identical  towing  forces.  Printout  of  cable  variables  is  desired  only  at 
end  points  of  cable  segments.  The  characteristics  of  the  fluids,  towline,  and 
depressors  are  listed  in  Table  7,  program  output  is  listed  in  Table  8. 
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Table  3.  Physical  characteristics  of  sample  problem  1. 


Speed  *  V  (kn) 

10.0 

Fluid  Density,  p  ( lb-sec^/f t^) 

1.9905 

Cable 

Diameter,  d  (in.) 

0.35 

Weight  in  Fluid,  W  (lb/ft) 

0.15 

Ribbon  Loading 

Normal  Drag  Coefficient,  Cr 

1.25 

Side  Force  Coefficient,  Cs 

0.10 

Cable  Length,  Smax  (ft) 

300.0 

Bare  (Pode)  Loading 

Pode  Frictional  Parameter,  f 

0.01 

Normal  Drag  Coefficient,  Cr 

1.90 

Side  Force  Coefficient,  C 

s 

0.10 

Total  Depth,  Zstop  (ft) 

300.0 

Buoy 

Drag  Force,  Fx  (lb) 

-500.0 

Side  Force,  (lb) 

0.0 

Down  Force,  F_  (lb) 
z 

-2000.0 

I 


Tabla  4.  Program  output  fo*  sample  problem  1» 


l  ■  _ 

|  SUBMARINE  TOWINS  A  SUBMERSED  BUOY  USING  A  STEEL  CIRCULAR  CABLE 


U-10.00  Knots 


Rho*1 .9905  Lb»SecA2/Ft  *4 


RIBBON  L0A0IN6 


fn-  0.4986-0.2499  C0S<Phi >+0.2527  SIN(Phi >-0.2487 
ft— 0.2255+0.3417  C0S(Phi  >+0.2255  SIN<  Phi  >+0.0000 
fs-  1.4792+0.0841  C0S(Phi >-3.0425  SIN< Phi  )-1 .5633 


C0S<2*Phi >+0.0000  SIN( 2*Phi  ) 
C0S(2«Phi >-0.081 1  SIN( 2*Phi  ) 
C0S( 2»Phi >+1 .5213  SIN(2*Phi) 


Cr-  1 .250 
Cs-  .100 
d-  .350  In. 

W-  .1S0  Lb/Ft 


BODY  FORCE  COMPONENTS 
Fx-  -500.00  Lb 
Fy-  0.00  Lb 
Fz-  -2000.00  Lb 


Sref ( Ft > 

Sstr( Ft  ) 

X<  Ft  > 

Y(  Ft  > 

Z(  Ft  > 

Phi ( Deo ) 

Beta( Deg  > 

T\  Lbs  > 

0.0 

0.0 

0.0 

0.0 

0.0 

76.0 

180.0 

2061.6 

50.0 

50.0 

17.2 

.4 

46.8 

63.9 

181  .0 

2081 . 1 

100.0 

100.0 

43.1 

1  .6 

89.5 

53.9 

182.3 

2112.2 

150.0 

150.0 

75.4 

3.6 

127.6 

45.8 

183.7 

2151  .6 

200.0 

200.0 

112.2- 

6.2 

161  .3 

39.4 

185.2 

2197.  1 

250.0 

250.0 

152.3 

9,3 

191  .0 

34.2 

186.5 

2247.2 

300.0 

300.0 

194.7 

12.6 

217.4 

30.0 

187.8 

2300.5 

Rho-1.9905  Lb*Sec'2/Ft  A4 
PODE  LOADING 


fn-  0.5000+0.0000  COS<Phi >+0,0000  S IN< Phi >-0.5000  C0S< 2*Phi >+0 . 0000  SIN(2*Phi> 
ft-  0.0100+0.0000  COS<  Phi )+0 . 0000  SINIPhi >+0.0000  C0S(2*Phi  >+0.0000  5IN<2*Phi> 
fs-  1.4792+0.0841  CQS< Ph i >-3 . 0425  SIN( Phi >- 1 . 5633  COS< 2*Phi >+ 1 . 52 1 3  SIN<2*Phi> 


Cr-  1 

.900 

Cs- 

.100 

d- 

.350  In. 

W- 

. 150  Lb/Ft 

Sref ( Ft > 

Sstr<  Ft  > 

X ( Ft  > 

Y(  Ft  > 

300.0 

300.0 

194.7 

12.6 

350.0 

350.0 

238.9 

16.0 

400.0 

400.0 

284.7 

19.3 

450 . 0 

450.0 

331.4 

22.6 

500.0 

500.0 

378.8 

25.7 

523.0 

523.0 

400.8 

27.2 

FINAL  FORCE  COMPONENTS 

Fx-  - 

2226.63  Lb 

Fy- 

-142.75  Lb 

Fz- 

-647,27  Lb 

Z  <  Ft  > 

Phi ( OeQ  ) 

Betel  Deo  > 

T( Lbs > 

2  1  7 . 4 

30.0 

187.8 

2200.5 

240.4 

25.5 

189.0 

2304.9 

260.4 

->n  7 

tea 

190.0 

2309.8 

277.9 

19.5 

191.1 

2315.0 

293.4 

17.4 

192.0 

2320.5 

300.0 

16.6 

192.4 

2323.2 

I 
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Tabic  5.  Physical  characteristics  of  sample  problem  2. 


Speed,  V  (kn) 

2S.0 

Water  Density,  pw  (ib-sec^/ft*) 

1.9905 

Air  Density,  p4  (lb-sec^/f t*) 

0.0024 

Towl ine 

Normal  Drag  Coefficient,  Cr 

0.15 

Side  Force  Coefficient,  Cs 

0.0 

Thickness,  d  (in.) 

0.5 

Weight  in  Water,  Ww  d(lb/ft) 

0.2 

We^jjht  in  Air,  W^  (lb/ft) 

0.55 

Modulus  of  Elas"icity,  E  (lb/in. ^) 

9.0  x  106 

A 

Structural  Cross-Sectional  Area,  Ac  (in.4) 

0.1 

Normal  Loading  Function,  fR 

-4.102  ♦  4.330  cost 
+4.878  sint  -  0.228  cos2t 
-2.165  sin2t 

Tangential  Loading  Function,  ft 

0.078  +  0.383  cost 
-0.078  sint 

Side  Loading  Function,  fg 

1.000  sint 

Athwartship  Distance,  Ystop  (ft) 

500.0 

Tow  Height,  Zstop  -  Depth  (ft) 

15.0 

Paravane 

Drag  Force,  Fx  (lb) 

-1000.0 

Side  Force,  Fy  (lb) 

5000.0 

Down  Force,  F  (lb) 
z 

1000.0 

30 


Table  6.  Program  output  for  sample  problem  2 


SHIP  TOWING  A  PARAVANE  USING  A  FAIRED  FIBERGLASS  TOWLINE 


l'*25 . 00  Knots 
Rho-1.9305  Lb*SecA2/FtA4 

fn— 4. 1020+4.3300  C0S( Phi  )+4 . 8780  SIN' Phi >-0. 2280  C0S< 2*Phi >-2 . 1650  SIN<2*Phi> 
ft-  0.0780+0.3830  C0S( Phi  )-0 . 0780  SIN< Phi >+0 . 0000  COSi 2*Phi  )+0.0000  SIN(2*phi) 
fa-  0.0000+0.0000  COSIPhi  )+l .0000  SIN< Phi >+0 . 0000  C0S( 2*Phi  )+0.0000  SIN(2*Phi> 

Cr-  .150 
Cs-  0.000 
d-  .500  In. 

W-  .200  Kb /Ft 
E-  9.000  *  1 0A6  Lb/In. A2 
Ac-  .100  In . "2 

BODY  FORCE  COMPONENTS 
Fx-  -1000.00  Lb 
Fy-  5000.00  Lb 
F?«  1000.00  Lb 


Sref  <  Ft ) 

SstrtFt > 

X(  Ft  ) 

Y<  Ft  ) 

Z(  Ft  ) 

Phi ( Deg  ) 

Beta<  Deg  ) 

T( Lbs  ) 

0.0 

0.0 

0.0 

0.0 

0.0 

78.9 

78.7 

5196.2 

50.0  . 

50.3 

12.0 

-47.9 

-9.6 

73.6 

78.6 

5249.9 

100.0 

100.6 

26.1 

-94.6 

-19.  1 

69.0 

78.5 

5322.3 

150.0 

150.9 

4‘<  .8 

-139.9 

-28.4 

65.0 

78.4 

5410.9 

200.0 

201  .2 

70.4 

-183.9 

-37.5 

61.5 

78.2 

5513.2 

250.0 

251.5 

95.6 

-226.5 

-46.4 

58.5 

73.1 

5627.3 

300 . 0 

301.8 

122.9 

-267.9 

-55.2 

55.9 

78.0 

5751.2 

350.0 

352.1 

151  .9 

-308. 1 

-63.8 

53.6 

77.9 

5883.7 

400 . 0 

402.5 

182.5 

-347.2 

-72.2 

51.6 

77.8 

6023.4 

450.0 

452.8 

214.3 

-385 . 3 

-80.5 

49.9 

77.7 

6169.5 

500.0 

503.2 

247.3 

-422.4 

-88.7 

48.3 

77.5 

6321.1 

515.5 

518.8 

257.7 

-433.8 

-91  .2 

47.9 

77.5 

6369.2 

Rho-0.0024  Lb*SecA“/Ft A4 

fn— 4. 1020+4.3300  C0S( Phi  >  +  4 . 8780  S IN< Phi  )-0 . 2280  C0S(  2*Phi >-2 . 1 650  SIN(2*Phi> 
ft-  0.0780+0.3830  COS<  Phi  )-0 . 0780  SIN( Phi  )+0. 0000  C0S( 2< Phi  )+0 . 0000  SIN<2*Phi> 
fs-  0.0000+0.0000  COS(Phi  )+1 .0000  SIN(Phi >+0.0000  C0S( 2*Phi  >+0.0000  SIN(2*Phi> 


Cr=  .150 
Cs-  0.000 
d=  .500  In. 


W- 

.550  Lb /Ft 

E-  9 

.000  ♦  1  0A6 

Lb /In. "2 

Ac- 

.100  In. A2 

Sref < Ft  ) 

Sstr( Ft ) 

X  ( Ft  ) 

Y(  Ft  ) 

515.5 

518.8 

257.7 

-433 . 3 

550.0 

553.5 

281.0 

-458.9 

600.0 

603.9 

314.8 

-495.3 

606.4 

610.3 

3.9.1 

-500.0 

FINAL  FORCE  COMPONENTS 
Fx-  -4274.22  Lb 
Fy-  4610.46  Lb 
Fz«  1072.54  Lb 


Z<  Ft  ) 

Phi ( Deg ) 

Beta( Deg  ) 

T(  Lbs  ) 

-91.2 

47.9 

77.5 

6369.2 

-96.8 

47.9 

77.3 

6372.4 

-105.2 

47.9 

75.9 

6377. 1 

-106.2 

47.9 

76 . 9 

6377.7 

Table  7.  Physical  characteristics  of  sample  problem  3. 


Speed,  V  (kn) 

20.0 

Water  Density,  pw  (Ib-sec^/ft^) 

1.9905 

Air  Density,  oa  (lb-sec^/ft^) 

0 .0024 

Towline  (Pode  Loading) 

Pode  Frictional  Parameter,  f 

0.007 

Normal  Drag  Coefficient,  Cr 

1.80 

Side  Force  Coefficient,  Cg 

0.0 

Diameter,  d  (in.) 

0.625 

Weight  in  Water,  Ww  d(lb/ft) 

0.000 

Weight  in  Air,  Wfl  (Ib/ft) 

0.130 

Modulus  of  Elasticity,  E  (lb/in.^) 

23,367.0  +  14.4  (T/Ac) 

Cross-Sectional  Area,  Ac  (in.  ) 

0.307 

Distance  Between  Depressors,  Smaxl  (ft) 

50.0 

Depth  to  Lower  Depressor,  Zstopl  (ft) 

75.0 

Total  Vertical  Separation,  Zstop2  (ft) 

175.0 

Paravane 

Drag  Force,  Fx  (lb) 

-300.0 

Side  Force,  F^  (lb) 

0.0 

Down  Force,  F  (lb) 

Ck 

1500.0 

Table  8.  Program  output  for  sample  problem  3 


HELICOPTER  TOWING  TWO  DEPRESSORS  USING  A  BRAIDED  NYLON  TOWLINE 


0-20. 00  Knots 
Rho-1.9905  Lb*SecA2/Ft  '4 
PODE  LOADING 

fn-  0.5000+0.0000  C0S<  Phi )+0 . 0000  SIN< Phi >-0.5000  COSL 2*Phi >+0.0000  SIN(2*Phi> 
ft-  0.0070+0.0000  C0S( Phi >+0.0000  SINLPhi >+0.0000  COSL 2*Phi  )+0 . 0000  SIN(2*Phi> 
fs-  1.4792+0.0841  COSL  Phi  )-3 .0425  SIN<  Phi  >-1  .5633  C0S(  2*Phi  >+1  .521  3  SINL2*Phi> 

Or*  1 . 800 
Cs*  0.000 
d-  .625  In. 

W-  0.000  Lb/Ft 

E-L  23.3870  +.01 44*T/Ac  >  *  1 0 ' 3  Lb/In. A2 

Ac-  .307  In. *2 

BODY  FORCE  COMPONENTS 
Fx«  -300.00  Lb 
Fy»  0.00  Lb 
Fz-  1500.00  Lb 


Sref ( Ft  ) 

SstrL Ft  ) 

X(  Ft  > 

Y<  Ft  > 

Z(  Ft  > 

Phi ( Deg > 

BetaL Deg  ) 

T( Lb  s  > 

0.0 

0.0 

0.0 

0.0 

0.0 

78.7 

0.0 

1529.7 

50.0 

54.9 

44.4 

0.0 

-28.1 

14.7 

0.0 

1568.7 

BODY  FORCE  COMPONENTS 

Fx« 

-300.00  Lb 

Fy* 

0.00  Lb 

Fz- 

1500.00  Lb 

Sref< Ft  ) 

Sstr<  Ft  > 

X  v  r  t  ) 

Y<  Ft  > 

Z<  Ft  > 

Phi ( Deg  > 

Beta( Deg  > 

T( Lbs > 

50.0 

54.9 

44.4 

0.0 

-28.1 

46.2 

0.0 

2627.5 

185.6 

209.1 

1R9.3 

0.0 

-75.0 

8.4 

0.0 

2735.6 

Rho-0 . 0024  Lb*SecA2/Ft  "4 


PODE  LOADING 

fn-  0.5000+0.0000  COSL Phi >+0.0000  SINLPhi >-0.5000  COSL 2*Phi >+0 . 0000  SIN(2*Phi> 
ft-  0.0070+0.0000  COS< Phi  >+0.0000  SIM  Phi  >+0.0000  COSL 2*Phi >+0.0000  SINL2*Phi> 
fs-  1.4792+0.0341  COSLPhi >-3.0425  SIN( Phi  >- 1 . 5633  COS ( 2*Phi >+ 1 . 52 1 3  SINL2*Phi> 


Cr»  1 . 800 
Cs*  0.000 
d“  .625  In. 

W-  .130  Lb/Ft 

E»<  23.3670  +  . 0 1 44*T /Ac  )  *10' 3  Lb/In. "2 

Ac-  .307  In.  A2 


Sref ( Ft  ) 

Sstr ( Ft  ) 

XL  Ft  ) 

Y<  Ft  > 

Z<  Ft  > 

Phi ( Deg  ) 

BetaL  Deg  > 

T( Lbs > 

185. 6 

209.1 

189.3 

0.0 

-75.0 

8.4 

0.0 

2735.6 

746.9 

842.2 

814.4 

0.0 

-175.0 

9,8 

0.0 

2747.6 

FINAL  FORCE  COMPONENTS 
Fx»  -2707.38  Lb 
Fy-  0.00  Lb 
Fz-  463.5!  Lb 
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R.  KNUTSON 


t 1/13/85 


10  I 
20  I 

30  I  Program  CABL36  calculates  the  three-dimensional 

40  I  CA8L36  static  configuration  of  an  extensible,  flexible 

50  I  cable  in  a  uniform  stream, 

60  I 

70  Jr - 

80  ! 

90  I  Main  input  section 

1001 

1  101 - 

120  RAD 

130  DIM  At( 801 

140  111;  Output-0 
150  Continue-0 

160  S-0 

170  Ss-0 

180  X-0 

1 90  Y-0 

200  Z-0 

210  T-0 

220  Fx-0 

230  Fy-0 

240  Fz-0 

250  INPUT  'FOR  HARDCOPY,  ENTER:  1  IFor  CRT,  press  CONTI*, Output 

260  IF  Output-0  THEN  PRINTER  IS  1 

270  IF  Output-1  THEN  PRINTER  IS  701 

280  INPUT  “ENTER:  TITLE* ,A» 

290  PRINT  * - — 


300  PRINT  A* 

310  PRINT  “- 


320  PRINT 

330  INPUT  “ENTER  SPEED:  U  ( Knots  )“,U 

340  PRINT  USING  1 1 1 "V-“ ,U , “Knots" 

350  II :  IMAGE  2A.DD.20 .X  ,5A 
360  M2:  PRINT 

370  INPUT  “ENTER  FLUID  MASS  DENSITY:  Rho  (Lb-Sec^/Ft^  )'  ,Rho 

3e0  PRINT  USING  I2i“Rho-“ ,Rho . *Lb«Sec*2/Ft‘4“ 

390  12:  IMA6E  4A ,Z . 40 ,X . 1 3A 

400  INPUT  “ENTER  CABLE  LOADING:  Pode-ti  Ribbon-2l  Other-3*  .Loading 

410  IF  Loading-1  THEN  50SUB  Pode 

420  IF  Loadlng-2  THEN  GOSUB  Ribbon 

430  IF  Loading-3  THEN  GOSUB  Other_loadi ng 

440  PRINT 

450  PRINT  USING  13 1 “ fnu" ,Fa0 ,Fal , "COS (Phi  )*  ,Fb I  , *SIN<  Phi >“  ,Fa2  ,  “C0S( 2 -Phi  )“ 

,Fb2,“SIN(2*Phi  )“ 


460  PRINT  USING  13 1”  ft-"  ,Ga0  ,Ga  I  .  “COS(  Phi  )“  ,Gb  1  ,  ”SIN<  Phi  >  ’’ 

,Gb2  ,  “SIN<  2»Phi  >“ 

470  PRINT  USING  13 i “ f s-“ ,Ha0 ,Ha 1 , “COSiPhi >* ,Hb 1 , * SIN< Phi > “ 

,Hb2  ,  “SIN!  2 'Phi  )“ 


480  13: 
490 
500 
510 

520  14: 

530 

540 

550  15: 


IMAGE  3A  ,MZ  .40  ,2(  SZ .  4D  ,X  ,3A  )  ,2<  SZ .  4D  ,X  ,  10A  ) 
PRINT  ' 

INPUT  “ENTER  CABLE  NORMAL  CRAG  COEFF:  Cr*  ,Cr 
PRINT  USING  14 i “Cr-“ ,Cr 
IMAGE  4X ,3A  ,00 . 3D 

INPUT  “ENTER  CABLE  SIOE  FORCE  COEFF:  Cs“,Cs 
PRINT  USING  I5t  “Cs-“  ,Cs 
IMA6E  4X  ,3A  ,M0 . 30 


,Ga2  ,"C0S(2*Phi )“ 
,Ha2  ,"C0S(2*Phi  )“ 
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560  INPUT  "ENTER  CABLE  THICKNESS:  d  (In.)". Ola 

570  PRINT  USIN6  16 i" d«" .Oia . “ In . “ 

580  16:  IMAGE  4X  ,2A  ,30 .30  ,X  ,3A 

590  INPUT  "ENTER  CABLE  WEIGHT  IN  FLUID:  W  (Lb/Ft  )".W 

600  PRINT  USING  17 » “W-" ,W."Lb/Ft“ 

610  17:  IMAGE  4X  ,2A  ,MDQ. 3D ,X  ,5A 

620  INPUT  "ENTER  CABLE  ELASTICITY:  Non-Elast ic-1 t  Elast lc-2 " .Elas 

630  IF  Elas<>2  THEN  M3 

640  INPUT  "ENTER  MODULUS  OF  ELAS:  E-E0+EI»T/Ac  ( Lb/In. "2 >“  ,E0 ,E1 

650  IF  EI-0  THEN  PRINT  USIN6  18 1 "E-*  ,E0/ 1 000000 . " » 1 0*6  Lb/In. *2" 

660  18:  IMAGE  4X  ,2A  ,30 . 30  ,X  ,  UA 

670  IF  E1O0  THEN  PRINT  USING  IB  i  "E-(  *  ,E0/ 1 000  ,E1  / 1  000  , "  *T/Ac  >  *  1 0  ~3  LG /In 

.  "2  " 

680  19:  IMAGE  4X ,3A ,30 . 40 ,SDO .4D  ,2 1 A 

690  INPUT  "ENTER  CABLE  CROSS-SECTIONAL  AREA:  Ac  (In."2)*,Ac 

700  PRINT  USING  1 1 0t “Ac-" .Ac , “ In. A2“ 

710  110:  IMAGE  4X ,3A ,00 . 30 ,X ,5A 
720  M3:  IF  Continue-2  THEN  M5 

730  M4 :  INPUT  "ENTER  BOOY  FORCE  COMPONENTS:  Fx  ,Fy  ,Fz  (Lb  >“  ,Fxo  ,Fyo  ,Fzo 
740  PRINT 

750  PRINT  "BOOY  FORCE  COMPONENTS" 

760  PRINT  USING  1 1  I | "Fx-" .Fro ."Lb" 

770  PRINT  USING  1 1  1 : "Fy-" ,Fyo  .  "Lb” 

780  PRINT  USING  1 1  I t "Fz-“ .Fzo , "Lb" 

790  III:  IMAGE  4X  ,3A  ,MSO . 20 .X  ,2A 

000  MS:  INPUT  "ENTER  MAXIMUM  CABLE  LEN6TH:  Snax  (Ft)",Smax 

810  INPUT  "ENTER  CABLE  LENGTH  PRINTOUT  INTERVAL:  Sine  (Ft)", Sine 

820  Cstop-0 

830  Xstop- 1 0AS 

840  Ystop- I 0"5 

850  Zstop- 1 0*5 

860  Phidstop--I0"3 

870  INPUT  "ENTER  LIMIT:  Xstop  (Ft)  C  I  f  none,  press  CONTI"  ,Xsiop 

880  INPUT  "ENTER  LIMIT:  Ystop  (Ft)  [If  nono  ,  press  CONTI", Ystop 

890  INPUT  "ENTER  LIMIT:  Zstop  (Ft)  ['If  none,  press  CONTI", Zstop 

900  INPUT  “ENTER  LIMIT:  Phistop  ( Oeg )  ilf  none,  press  CONT ] “ .Phidst op 

910  P)u  s top-Phidst op/S7 . 2958 

920  PRINT 

930  PRINT  "  Sref(Ft)  Sstr(Ft)  X  <  Ft  )  Y<  Ft  )  Z  <  Ft  )  PhilDeg'  BetalDe 

g  )  T< Lbs  ) ” 

940  IF  Continued  THEN  G0SU8  Tangle 

950  GOSUB  Integrator 

960  Conti nue-0 

970  INPUT  "ADD  BOOY.  ENTER:  !|  ADD  CABLE,  ENTER:  2  [NEW  CASE,  press  CONTI", 

Cont 1 nue 

980  IF  Continue-I  THEN  M4 

990  IF  Continue-:  THEN  M2 

1000  PRINT 

1010  PRINT  “FINAL  FORCE  COMPONENTS" 

1020  PRINT  USING  1 1  I : “Fx-“ ,Fx  .  “Lb“ 

1030  PRINT  USING  111:  “Fy-"  ,Fy  ,"Lb“ 

1040  PRINT  USING  1 1  I  I  "Fz-“  ,Fz  ,  "Lb  " 

1050  IF  Output-1  THEN  PRINT  CHR*<12> 

1060  GOTO  Ml 
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1070 

1060 

1090 

1100 

1110 
1120 
1 130 
1 140 
1 150 
1 160 
1170 
1 180 
1190 
1200 
1210 
1220 
1230 
1240 
1750 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 


I - 

I 

Tangle:  I  Subroutine  ‘Tangle*  calculates  starting  cable 

I  angles  and  tension  from  input  force  conponents. 


Fx-Fx+Fxo 

Fy-Fy+Fyo 

Fz-F2+FzO 

Tzy-SQR(Fz»Fz+Fy«Fy ) 

T-SQR< Tzy»Tzy+Fx*Fx ) 

IF  Tzy-0  THEN  Beta-0 

IF  (Fz-0)  ANO  ( Fy<>0  )  THEN  Beta-PI *Fy/< 2*ABS< Fy  ) ) 

IF  Fz>0  THEN  Beta-ATN<  Fy/Fz  > 

IF  ( Fz<0 )  ANO  (Fy-0)  THEN  Beta-PI 

IF  ( Fz<0  )  AND  <  Fy<>0  )  THEN  Beta-ATN( Fy/Fz )+PI*Fy/ABS( Fy ) 
IF  <  Fx-0  )  AND  ( Tzy<>0  >  THEN  Phi-PI/2 
IF  Fx<0  THEN  Phi«ATN( -Tzy/Fx  ) 

IF  Fx>0  THEN  Phi-ATN< -Tzy/Fx >+PI 
IF  T-0  THEN  60SUB  Critang 
IF  T< . 1  THEN  T-.1 
RETURN 


I  Subroutine  “Critang"  calculates  critical  angles 

Critang:  1  to  begin  the  integration  when  the  starting 

I  tension  is  zero. 


I - 

Beta-0 

Phi-0 

IF  U-0  THEN  P3 
R-. 1 1 8723*Rho*Cr«Dia»U*U 
Fs». 1  1 8723»Rho«Cs»Dia»U*U 
Phlnc-.09 
PI:  Sp-SIN(Phi) 

Cp-COS( Phi  ) 

S2p»SIN( 2*Phi  ) 

C2p-C0S(2*Phi ) 

H-Fs*(Ha0+Ha1 *Cp+Hb 1 »Sp+HaZ*C2p+Hb2*SZp  > 

IF  ABS( U >>«ABS( H  )  THEN  PZ 
Phinc-Phinc/2 
Phi-Phi-Phinc 
GOIO  PI 

P2 :  IF  W>0  THEN  Beta=ASN( H/W  ) 

IF  (U<0)  AND  <  H>-0  )  THEN  Beta— ASN<  H/W  )-PI 
IF  <U<0>  ANO  <u/0>  7!--'  eta— ASN(H/U)+PI 

F— R* <  Fa0  -  M  -rbl'  2*C2p+Fb2»S2p  ) 

Sum“F+W»COS( Phi  )*COS\ Jeta  > 

IF  ABS< Sun  X. 00001  THEN  P3 
IF  Sum<0  THEN  Phinc=Phinc/2 
Phi-Phi+Phinc«Sun/ABS(  Sun ) 

60T0  PI 
P3 : RETURN 
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Subroutine  “Integrator"  solve*  the  differential 
aquations  of  cable  force  equilibrium  using  a 
self-starting,  modified  Euler's  method. 


1610  ! - 

1620  I 
1630  1 

1640  Integrator:  • 

1650  1 
1660  1 

1670  ! - 

IG80  Sinc«ABS(  Sine  ) 

1690  IF  Smax<0  THEN  Sinc--Sir.c 

1700  Dias-Dia 

1710  Ws-U 

1720  N I  :  Phid-Phi»57.2958 

1730  Betad-Beta*57.2958 

1740  PRINT  USING  1 1 2 1  S  ,Ss  ,X  ,Y  ,Z  ,Phid  .Betad  ,T 
1750  112: IMAGE  S< M5D . 0 ,X > ,2X ,2(  M3D .0  ,4X  )  ,7D . D 
1760  IF  (Cstop-I)  OR  <  ABS<  Smax-S  K-,01  )  THEN  N9 

1770  SI  "INT ( <  ABS<  S  >+ .  0 1  >/ABS< Sine  >  >*Sinc+Sinc 


1730  N7 : 
1790 
1800 

1810  N3 : 
1820 
1830 
1840 
1850 

1860  N4 : 
1870 
1880 
1890 
1900 
1910 
1920 
1930 
1940 
I960 
1960 
1970 
1980 

1990  N5 : 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
2110 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 
2200 


Dels-SI-S 

IF  ABS( De 1 s  )>50  THEN  Del s«50*Dels/ABS< Dels  ) 
IF  ABS(5+0els )>-ABS< Smax >  THEN  Dels-Smax-S 
Oelss-Oels 
1-0 

Delp-0 
Delb-0 
Oelt»0 
1  =  1  +  1 

Delp 1 -Dalp 
Delb  1 -Oalb 
Betav-Beta+0elb/2 
Phiav«Phi+0elp/2 
Tav-T+Delt/2 
IF  Elas< >2  THEN  N5 
E-E0+E1 *Tav/Ac 
IF  EI-0  THEN  B-1+Tav/(Ac*E > 

IF  E1O0  THEN  B-1+LOG<E/E0)/EI 
Oias»Dia*SQR< 1 /B ) 

Ws»W» (  1/B  ! 

Delss*Dels»8 

R- , 1 1 872j*Rho»Cr*0ias*V»V 
Fs-. 1  I 8723*Rho*Cs*0ias*V*V 
Sp-SIN<  Phiav  1 

IF  ABS(Sp  X.0001  THEN  Sp-.000l 
Cp-C0S<  Ph iav  ) 

S2p-S  IN(  2*Phiav  > 

Sb«SIN(Betav ) 

Cb*C0S( 8etav ) 

Phiav1-AB3(ASN(SIN< Phiav))  ) 

Spl-SIN(Phiav)  ) 

Cp I —COS < Ph lav  I  ) 

S2p  1  —  S IN <  2*Phiav  I  1 
C2p  1 -C0S(  2»Pruav  1  > 

FI--R»!Fs0+Fal»Cpt+Fb I *Spl+Fa2»C2pl+Fb2*52pl  ) 
G1— R«( Ga0+Ga I *001+661 *Spl +Ga2*C2p1+Gb2*S2p1  > 
HI -Fs*( Ha0+Ha I *Cp I +Hb I *Sp 1 +Ha2*C2p 1 +Hb2»52p I  > 
IF  Sp>-0  THEN  F-F I 
IF  Sp'0  THEN  F— FI 
IF  Cp  >-0  THEN  G-G1 
IF  Cp10  THEN  C»--G  1 

IF  (Loadmg<»2)  ANO  (S2p.>-0)  THEN  H-HI 
IF  ( Loading<-2  )  AND  <S2p<0>  THEN  H--HI 
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2210  IP  <  Loading>2  )  AND  <Sp>-0>  THEN  H-H1 

2220  IF  ( Loading>2  )  AND  <Sp<0>  THEN  H— HI 

2230  Delt-< -6+Ws»Sp»Cb >»Oalaa 

2240  Delb«< H~Wa*Sb  >*0el8a/<  Tav«Sp ) 

22S0  0elp-< F+Wa*Cp»Cb  )»0elsa/Tav 

2260  Delst-Oels 

2270  IF  I>10  THEN  Oela-Dalal/2 

2280  IF  A8S<0alb )>.09  THEN  Oel s-ABS< ;0?/Oelb XOels 1 

2290  IF  ABS< Delp )> .09  THEN  Oels-ABS< .07/Delp >*Oela1 

2300  IF  DelaODels  I  THEN  N3 

2310  IF  <  A8S< Delp-Dalpl )> .0002  >  OR  ( ABS(Oelb-Delb) )>.0002 >  THEN  N4 

2320  S^S+Dels 

2330  Sa-Sa+Oelas 

2340  X«X+Oalaa»Cp 

2350  Y-Y-0elas*Sp*Sb 

2360  Z-Z-Delaa*Sp*Cb 

2370  T-T+Dalt 

2380  Phi-Phi+Delp 

2390  Beta-Beta+Delb 

2400  *  IF  Catop-I  THEN  N7 

2410  IF  ABS(  X  )>ABS<  Xatop  )  THEN  Nt» 

2420  IF  ABS( Y  »ABS( Ystop )  THEN  N6 

2430  IF  ABS<  Z  )>ABS<  Zatop  )  THEN  N6 

2440  IF  Phi<Phiatop  THEN  N6 

2450  IF  < ABS< Snax-S  )> .01  )  AND  < ABS< S 1 -S » . 01  >  THEN  N2 

2460  SOTO  N I 

2470  NS:  Catop-l 

2480  N7 :  IF  ABS<  ABS<  X  >--A8S(  Xatop  )  X  .01  THEN  N1 

2490  IF  ABS<  A8S<  Y  )-ABS( Yatop ) >< .01  THEN  N1 

2500  IF  ABS( ABS<  Z )-ABS( Zstop  >  X.01  THEN  N1 

2510  IF  ABS(Phi-Phiatop K.0005  THEi.  N1 

2520  Oeis"ABS<  Dels ) 

2530  IF  ABS<X  )>ABS( Xatop )  THEN  N8 

2540  IF  ABS( Y)>ABS( Yatop >  THEN  N8 

2550  IF  ABS<  Z  )>ABS( Zatop  )  THEN  N8 

2560  IF  Phi<Phiatop  THEN  N8 

2570  Dela-Oels/2 

2580  IF  S<0  THEN  Dels— Dels 

2590  GOTO  N3 

2600  N8 :  Dels— Dela/2 

2610  IF  S<0  THEN  Dels— Dels 

2620  GOTO  N3 

2630  N9:  Fx— T*C0S(Phi  > 

2640  Fy"T*SIN(Phj.  )*SIN(Beta) 

2650  Fz-T*SIN(Phi )*C0S<6eta ) 

2660  RETURN 
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2670  ! 
2660  ! 


2690 

2703 

2710 

2720 

Pod«:  1 

I 

\ 

i - - - 

Subroutine  “Pcda"  providao  coefficients  for 
Poda  or  bara  cable  loading  function*. 

2730 

Fa0-,5 

2740 

Fa  l»0 

2750 

Fbt-0 

2760 

Fa2— .5 

2770 

Fb2-0  ' 

2760 

6a  1-0 

2",90 

6b  1  -0 

2800 

6«2-0 

2810 

6b2-0 

2820 

Ha0-l .4792 

2830 

Hal-. 0841 

2840 

Hb  1—3.0425 

2850 

Ha2— 1 .5633 

2860 

Hb2- 1 .5213 

2870 

PRINT 

2880 

PRINT  "PODE 

LOADING” 

2890 

INPUT  “ENTER 

PODE  FRICTIONAL  PARAMETER:  f“.ba0 

2900 

RETURN 

2920 

( 

2930 

Ribbon:  1 

Subroutine  “Ribbon"  provides  coefficients  for 

2940 

! 

ribbon-faired  cable  loading  functions. 

2950 

1 

£960 

2970 

Fa0- . 4986 

2980 

Fal— .2499 

2990 

Fbl-,2527 

3000 

Fa2— .2487 

3010 

Fo2*0 

3020 

6a0— .2255 

3030 

Gal-. 3417 

3040 

6b I-. 2255 

3050 

6a2-0 

3060 

6b2—  .081 1 

307C 

Ha0»1 .4792 

3080 

Hal -.0841 

3090 

Hbl—  3.0425 

3100 

Ha2— 1 .5633 

3110 

Hb2-1 .5213 

3120 

PRINT 

3130 

PRINT  “RIBBON  L0ADIN6" 

7140 

RETURN 

XLS0 

3160 

1 

3170 

Other_loading:  • 

Subroutine  "Other _ loading"  allows  input  of 

3180 

1 

generalised  cable  loading  functions. 

3190 

1 

3200 

3210 

INPUT  “ENTER 

NORMAL  L0A0IN6  FUNC  .  fn:  A0  ,A 1  ,B I  ,A2  ,B2 "  .Fa0  .Fa  1  .Fb 1  Fa2  F 

b2 

3220 

INPUT  “ENTER 

TANGENTIAL 

LOADING  FUNC,  ft:  A0  ,  At  ,B  1  ,  A2  ,62 '  ,Ga0,Gal  6b  I  6 

a2  ,Gb 

2 

3230 

INPUT  “ENTER 

SIDE  LOADING  FUNC,  fs:  A0  ,  A 1  ,B  1  ,  A2  .82  “  ,Ha0  ,Ha  1  .Hb  1  ,Ha2  .Hb2 

3240  RETURN 
3250  END 
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